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1. Introduction 

Trace formulae provide one of the most elegant descriptions of the classical-quantum corre- 
spondence. One side of a formula is given by a trace of a quantum object, typically derived 
from a quantum Hamiltonian, and the other side is described in terms of closed orbits of the 
corresponding classical Hamiltonian. In algebraic situations, such as the original Selberg trace 
formula, the identities are exact, while in general they hold only in semi-classical or high-energy 
limits. We refer to a recent survey |l4| ] for an introduction and references. 

In this paper we present an intermediate trace formula in which the original trace is expressed 
in terms of traces of quantum monodromy operators directly related to the classical dynamics. 
The usual trace formulae follow and in addition this approach allows handling effective Hamilto- 
nians. 

Let P = (l/i)hd x be the semi-classical differentiation operator on the circle, x G S 1 = R/27rZ, 
< h < 1. The classical Poisson formula can be written as follows: if / € C£°(R) then 

(1.1) tr/(PM) = JL ^ Jf {z/h) (e^) k ± ( e >W*) dz , 

\k\<N JK 

where N depends on the support of /, and we think of M(z, h) — e 27Tlz ^ h : C — > C as the 
monodromy operator for the solutions of P — z. It acts on functions in one dimension lower (zero 
dimension here) , identified geometrically with the functions on the transversal to the closed curve 
(S 1 here), and analytically with ker (P — z) (C here). 

Now let P be a semi-classical, self-andjoint, principal type operator, with symbol p (for instance 
P = — h 2 A + V(x), p = £ 2 + V(x)), and let 7 C p _1 (0) be a closed primitive orbit of the Hamilton 
flow of p. We can define the monodromy operator, M(z, h) for P — z along 7, acting on functions 
in one dimension lower, that is, on functions on the transversal to 7 in the base. We then have 

Theorem 1. Suppose that there exists a neighbourhood of"/, fl, satisfying the condition 

(1.2) men and exptH p (m) = m , p{m) = , < \t\ < TN => m E 7 , 

where T is the primitive period of 7. If f G C C °°(R), supp / C (-NT + C, NT - C) \ {0}, 
C = C(p) > 0, x G C^°(R) ; and A G (X) is a microlocal cut-off to a sufficiently small 
neighbourhood 0/7, then 

1 r d 

(1.3) tr f(P/h) X (P)A=— ]T tr f{z/h)M{z,h) k —M{z,h)x(z)dz + 0(h°°), 

111 -JV-1 R z 
where M(z, h) is the semi- classical monodromy operator associated to 7. 

The dynamical assumption on the operator means that in a neighbourhood of 7 there are no 
other closed orbits of period less than TN, on the energy surface p = 0. We avoid a neighbourhood 
of in the support of / to avoid the dependence on the microlocal cut-off A. 
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The monodromy operator quantizes the Poincare map for 7 and its geometric analysis gives 
the now standard trace formulae of Selberg, Gutzwiller and Duistermaat-Guillemin (see [fy for a 
recent proof and a historical discussion, and Sect.]?] for a derivation based on Theorem [l]). The 
term k = — 1 corresponds to the contributions from "not moving at all" and the other terms to 
contributions from going \k + 1| times around 7, in the positive direction when k > 0, and in the 
negative direction, when k < — 1. For non-degenerate orbits we analyse the traces on monodromy 
operators in Sect.^ and recover the usual semi-classical trace formulas in our general setting - see 
Theorem |^. 

Theorem |l| is a special case of the more general Theorem || presented in Sect|| Motivated 
by effective Hamiltonians in which the spectral parameter appears non-linearly, we give there a 
trace formula for a family P(z) with the special case corresponding to P — z. For an example of 
a use of effective Hamiltonians in an interesting physical situation we refer to Q|. The effective 
Hamiltonian described there comes from the "Peierls substitution" , and the celebrated "Onsager 
rule" is a consequence of a calculation of traces. 

The point of view taken here is purely semi-classical but when translated to the special case 
of C°°-singularities/high energy regime, it is close to that of Marvizi-Melrose Jl0| and Popov |T^| . 
In those works the trace of the wave group was reduced to the study of a trace of an operator 
quantizing the Poincare map. In jl2| it was used to determine contributions of degenerate orbits 
and our formula could be used for that as well. 

Acknowledgments. We would like to thank the Erwin Schrodinger Institute, and the France- 
Berkeley Fund for providing partial support for this work. The second author gratefully acknowl- 
edges the partial support of the National Science Foundation under the grant DMS-9970614. 



2. Outline of the proof 

T o present the idea of our proof we use it to derive the classical Poisson summation formula 
(1.1). The left hand side there can be written using the usual functional calculus based on 
Cauchy's formula: 




= R±iR, 



where we take the positive orientation of R and R > is an arbitrary constant. We make an 
assumption on the support of the Fourier transform on /: 

(2.2) supp / C (-2nN, 2ttN) . 

We would like to replace (P — z)^ 1 by an effective Hamiltonian which measures the obstruction 
to the solvability of (P — z)u = f. For that we introduce a Grushin problem (see for instance || 
for applications of this method in spectral problems, and for references): 

(2.3) V{z) ^ ( ) : H\&) x C -> ^(S 1 ) x C, 

where R±(z) should be chosen so that V{z) is invertible. If we put 

R+u = u(0) , 



then we can locally solve 



f (P-z)u = 
\ R+u = v 
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by putting 

u = I+(z)v = exp(izx/h)v , — e < x < 2ir — 2e . 
This is the forward solution, and we can also define the backward one by 

u = I-(z)v = exp(izx/h)v , — 2ir + 2e < xe . 
The monodromy operator M(z, h) : C — ► C, can be defined by 
(2.4) I + (z)v(n)=I-(z)M(z,h)v('iv), 
and we immediately see that 

M(z,h) = exp— — . 

We use I±(z) and the point ir to work with objects defined on S 1 rather than on its cover: a 
more intuitive definition of M(z, h) can be given by looking at a value of the solution after going 
around the circle but that has some technical disadvantages. 
Let x G [0, 1]) have the properties 

X(x) = 1 , -e<x<n + e, x(ir)=0, —it + 2e < x < -It , 

and put 

E + (z) = X I + (z) + (1 - x)I-(z) ■ 

We see that 

(P - z)E + = [P, X }I+(z) - [P, X ]I-(z) - [P,x]-I+{z) - [P,x]-I-(z) , 

where [P, x]— denotes the part of the commutator supported near tt. This can be simplified using 
©: 

(P - z)E+ + [P, X ]-I-(z)(I - M(z, h)) = 0, 

which suggests putting 

R-(z) = [P lX ]-I-(z) 1 



so that the problem 



has a solution: 



(P - z)u + R-(z)u- 
R+(z)u = v 



u = E + (z)v 
U-=E- + {z)v ' 



with P_ + (z) = I -M(z,h). 

One can showQ that with this choice of R±(z), ( |2.3| ) is invertible and then 



P(z)- 1 = £(z)=( E W E+{Z) ) 
I Kz) y E _^ E_ + {z) ) 



where all the entries are holomorphic in z, and E+(z), E \-(z), are as above. The operator 

E \-{z) is the effective Hamiltonian in the sense that its invertibility controls the existence of 

the resolvent: 

(2.5) (P - z)- 1 = E{z) - E + (z)E-+(z)- 1 E-(z) . 



1 In this situation it is quite easy but it will be done in greater generality in Sect.^. 



4 



J. SJOSTRAND AND M. ZWORSKI 



Inserting this in (2.1) and using the holomorphy of E(z) gives 





=--/ 


(x) 


2m J T 




= --/ 




2m J T 



where we used the cyclicity of the trace. Differentiating £(z)V(z) = Id shows 

E_{z)E+{z) = d z E_ + (z) + E_(z)d z R-(z)E_+(z) , 
which inserted in the previous identity gives 



where we eliminated the other term using countour deformation. 
We now use the expression for E |_ to write 



P 



1 



/ (£) tr d z M{z, h)(I - M(z, ft)) -1 dz 



+ — — r / f(^-)trd z M(z,h)M(z,h)- 1 (I-M(z,hy 1 y 1 dz. 



The assumption (2.2) and the Paley- Wiener theorem give 



\f(z/h)\ < e 27rJv l Im2 l/' 1 (Rez/^)- oc \ 

Writing 

JV-l 

(/ - M(z, h))- 1 = M ( z > h ) k + M ( z > ~ M ( z > 1 

for r + , and 

JV 

M(z, h)-\l - M(z, h)- 1 )- 1 = Y M ( z > h V k + M ( z > hy 1 *- 1 ^ - M(z, hj)- 1 , 

fc=i 

for r_, we can eliminate the last terms by deforming the contours to imaginary infinities (R — > 
oo), and this gives (1.1). 

In the general situation we proceed similarly but now microlocally in a neighbourhood of the 
closed orbit described in Theorem [j] - see Sect. | for a precise definition of microlocalization. The 
formula (^[l]) has to be replaced by 

(2.6) tr / X (P)A = -!£/(£) B z x(z)(P - z^A C{dz) , 

where x is an almost analytic extension of x, that is an extension satisfying d z x{z) = 0{\ Iva.z\°°) 
- see Sect.^J, and we want to procceed with a similar reduction to the effective Hamiltonian given 
in terms of the monodromy operator. 

To construct the monodromy operator we fix two different points on 7, mo, toi (corresponding 
to and it in the example), and their disjoint neighbourhoods, W+ and W- respectively. We 
then consider local kernels of P — z near Too and m\ (that is, sets of disctributions satisfying 
(P — z)u — near TOj's), ker mj . (P — z), j = 0, 1, with elements microlocally defined in W±. and 
the forward and backward solutions: 

I±(z) : ker mo (P-z) — > ker mi (P - z) . 
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We then define the quantum monodromy operator, M.{z) by 

I-(z)M(z) = I+(z) , M(z) : ker mo {P-z) — » ker mo (P-z). 

The operator P is assumed to be self-adjoint with respect to some inner product (•,•), and we 
define the quantum flux norm on ker mo (P — z) as follows^: let x be a microlocal cut-off function, 
with basic properties of the function x m the example. Roughly speaking \ should supported 
near 7 and be equal to one near the part of 7 between W+ and W-. We denote by [P, x]w + the 
part of the commutator supported in W+ , and put 

(u,v) QF d = ([(h/i)P,x]w+u,v) , 11,1)6 ker mo (P - z) . 



It is easy to check that this norm is independent of the choice of % - see the proof of Lemma 4.4. 
This independence leads to the unitarity of M(z): 

(M(z)u,M(z)u)qf = (u,u) QF , u G ker TOo (P - z) . 

For practical reasons we identify ker mo (P—z) with I?'(]R" _1 ), microlocally near (0, 0), and choose 
the identification so that the corresponding monodromy map is unitary (microlocally near (0, 0) 
where (0,0) corresponds to the closed orbit intersecting a transversal identified with T*M™ _1 ). 
This gives 

M{z,h) : V'iR"- 1 ) — ► V'{R n - v ), 

microlocally defined near (0, 0) (see Sect.| for a precise definition of this notion) and unitary 
there. This is the operator appearing in Theorem |l| and it shares many properties with its simple 
version e~xjp(2iriz/h) appearing for S . 

As shown in the example of the Poisson formula, traces can be expressed in terms of traces 

of effective Hamiltonians (E |_(z) there). Hence in our final formula, we replace P — z by a 

more general operator P(z), for which we do not demand holomorphy z but only that P{z) is 
self-adjoint for z real and that it is an almost analytic family of operators. In Theorem |^ in 
Sect|] we will compute the trace of 

-itr J f(z/h) B z [x(z) d z P(z) P(z)- 1 ] A C(dz) , 

which for P(z) = P — z reduces to ( |2.6| ). 

The only prerequisite to reading the paper is the basic calculus of semi-classical pseudodiffer- 
ential operators (see ||). In SectJ^ we review various aspects of semi-classical microlocal analysis 
needed here. In Sect we define the quantum time and quantum monodromy. Then in SectJ^ we 
follow the procedure described for S 1 to solve a Grushin problem allowing us to represent P(z) _1 
near a closed orbit. That is applied in the proof of the trace formula in Sect.^, and in Sect.0 we 
derive the more standard trace formula in the case of a non-degenerate orbit. 

3. Semi-classical operators and their almost analytic extensions 
Let X be a compact C°° manifold. We introduce the usual class of semi-classical symbols on 

X: 

S m - k (T*X) = {aeC°°(T*X x (0,1]) : \d^a(x,^;h)\ < C a ,ph- m (Z) k -M} , 



2 See pj for an earlier mathematical development of this basic quantum mechanical idea. 
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and the class corresponding pseudodifferential operators, fy 1 ?' (X), with the quantization and 
symbol maps: 

o h : #™' fc (X) — > S m ^{T*X)/S m - 1 ' k - 1 {T*X), 
with both maps surjective, and the usual properties 

o h (AoB) = o h {A)(T h {B), 

-> ^ q m ' k (X) °A S m ' k {T*X)/S m - 1 > k -\T*X) 0, 

a short exact sequence, and 

a h o Op™ : S m - k {T*X) — ► S m ' k (T*X)/S m ~ 1 ' k ~ 1 (T*X) , 

the natural projection map. The class of operators and the quantization map are defined locally 
using the definition on K n : 

(3.1) Ofl(a)u(x) = J J a (Z±*, e^-y^l h u{y)dyd^ , 

and we refer to || or [|l3) for a detailed discussion. We remark only that unlike the invariantly 
defined symbol map, ah, the quantization map Op™ can be chosen in many different ways. 

In this paper we consider pseudo-differential operators as acting on half-densities and con- 
sequently the symbols will also be considered as half-densities - see ||, Sect. 18.1] for a general 
introduction and the Appendix to this paper for a semi-classical discussion. For notational sim- 
plicity we supress the half-density notation. The only result we will need here is that in Weyl 
quantization, the symbol is well defined up to terms of order 0(h 2 ) - see Appendix. 

For a € S m ' k (T*X) we define 

ess-supp ft a C T*X U S*X , S*X d = (T*X \ 0) /R+ , 
where the usual R_|_ action is given by multiplication on the fibers: (x,^) i— > (a;,i£), as 
ess-supp^ a = 

Z{(x,0eT*X : 3e>0d^d^a(x',O^O(h oc ), d(x,x') + |f - f| < e} 
UC{(£,£) 6T*I\0 : 3e>0 d«dPa(x'^') = 0(h O0 {£y co ), 

d(x,x') + i/\e\ + m\-e/\e\\<e}/R + 

For A e *™' fe (A), then define 

WFh(A) = ess-supp h a , A — Op™ (a) , 
noting that, as usual, the definition does not depend on the choice of Op™. For 
ueC°°((0,l] h ;V'(X)), 3 N , h~ N °u is bounded in V'{X), 

we define 

WF h (u) = C{(x,0 : 3Ae*° h >°(X)o h (A)(x,0^0, Au e h°°C°°((0, l] h ; C°°(X))} . 

When u is not necessarily smooth we can give a definition analogous to that of ess-supp^ a. Since 
in this note we will be concerned with a purely semi-classical theory and deal only with compact 
subsets of T* X this definition is sufficient for our purposes (for more general definitions of wave 
front set which include this usual semi-classical definition, see pdj|). 
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To discuss almost analytic continuation of semi-classical pseudodifferential operators let us 
first recall the scalar case. For / G C tx> (M), an almost analytic extension of / is / G C°°(C) such 
that locally uniformly 

BJ(z) = 0(\Imzn, fh=f- 

The almost analytic extensions were introduced by Hormander and are unique up to 0(\ lmz\°°) 
terms (see [|[ Sect. 8] and references given there). 
Suppose now that 

A(x)€C°°(R x ;*™' k (X)) 

is a smooth family of pseudodifferential operators. We can then find a(x) £ C°°(M. X ; S m ' k (T*X)) 
such that A(x) — Op™(a(x)). We then define the almost analytic extension of the family A{x) as 

A(*) = Op£(a(z)), 

where a(z) G C°°(C Z ; S m ' k (T*Xj) is an almost analytic extension of a(x). To justify this defini- 
tion we need the following easy 

Lemma 3.1. If a(x) G C°°(M. X ; S m ' k (T*X)) then there exists an almost analytic extension of a 
satisfying 

a{z) G C°°(C Z ; S m ' k (T*X)) , d^B z a(z)(x, f; h) = 0(\ Im z|°°(£) fc ~ l/31 ) . 

We will also need certain aspects of the theory of semi-classical Fourier Integral Operators. 
Rather than review the full theory we will consider a special class, to which the general calculus 
reduces in local situations. Thus let A(t) be a smooth family of pseudodifferential operators, 
A(t) = Op£(a(*)), a(t) G C~([-l,l] t ;S°- M (T*X)), such that for all t, WF(A{t)) <e T*X. We 
then define a family of operators 

U(t) : L 2 {X)^L 2 (X), 

(3 2) 

hD t U(t) + U(t)A(t) =0, U(0) = U £^°(X). 

This is an example of a family of h- Fourier Integral Operators, U(t), associated to canonical 
transformations n(t), generated by the Hamilton vector fields H ao r t ), where the real valued ao(t) 
is the /i-principal symbol of A(t), 

j t K(t)(x,0 = (K(t))*(H aoit) ( X> 0), K (0)(x,O = (x,O, (x,OeT*X. 

All that we will need in this note is the Egorov theorem which can be proved directly from this 
definition: when U in |T|) is elliptic (that is \a{U )\ > c> on T*X, then for B G *™' fe (X) 

a{V{t)BU{t)) = («(t)) V(fl) , 

' /> "" > ' V{t)U{t) - I , U{t)V{t) -Ie ^-°°'- 00 (T*X) 

where the approximate inverse is constructed by taking 

hD t V(t)-A(t)V{t) = 0, V(0) = V Q , V Q U - 1, U V - 1 G -°°(T*X) , 

the existence of Vq being guaranteed by the ellipticity of Uq. The proof of ( |3.3| ) follows from 
writing B(t) — V(t)BU(t), so that, in view of the properties of V(t), 

hD t B{t) = [A(t),B{t)} mod f^ 00 , B(Q) = B . 

Since the symbol of the commutator is given by (h/i)H ao ^a(B(t)), ( |3.3| ) follows directly from 
the definition of n(t). 
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If U = U(l), say, and the graph of k(1) is denoted by C, we conform to the usual notation 
and write 

U€l° h (XxX;C), C = {(x,£;y,-Ti) : (x, £) = K(y, V )} , 

which means that U is an fe-Fourier Integral Operator associated to the canonical graphs C. 
Locally all h- Fourier Integral Operators associated to canonical graphs are of the form U(l) 
thanks to the following well known 

Lemma 3.2. Suppose that U\,TJ% are open neighbourhoods of (0,0) E T*W n , and k : U\ — > U2 
is a canonical transformation satisfying k(0,0) = (0,0). Then there exists a smooth family of 
canonical transformations Kt ■ U\ — » XJ%, < t < 1, satisfying kq — id, k\ = k, Kt(0,0) = (0,0). 

Proof. Since the symplectic group, Sp(n, R), is connected we can first deform k so that c?k(0, 0) = 
Id. Hence, near (0,0), ((x(K(y, 77)), £(k(v, rf)); y, rj) 1— > (x, rj) is surjective, and on the graph of 
k, y and 77 can be regarded as functions of x and rj. Since the symplectic forms, —d((y,dr])) 
d((£, dx)) are equal, their difference can be written locally as a differential: 

(y,dri) + (^dx) (b^^x^), #(0,0) = 0, 

so that k : ((f>L(x, 77), 77) 1— » (x , 4>' x (x , rj)) . We could now take as our family 

K t : (t<j>! n (x,T}) + (l-t)x > r))^(x ) t<f>x(x,Tj) + (l-t)ri). 

The two steps can be connected smoothly by making the deformations flat at their junction. □ 



The almost analytic continuation of a family of h- Fourier Integral Operators defined by (3.2) 
is obtained by means of the following 

Lemma 3.3. Suppose that U(t) is defined by ( |3.2j) a nd that A(z) is an almost analytic con- 
tinuation of the family A{t), as given by Lemma 3.1. Let U{z) = U(t + is) be the solution 
of 

(3.4) -hD s U(t + is) + U(t + is)A(t + is) = 0, U(t + is) \ s = = U(t) . 
1 

Then for |Imz < h\ogh~ L we have 

(3.5) ||C/(z)|j i 2^ L 2 < Cexp(C|Imz|//7) 

(3.6) \\d z U(z)\\ L ^ L2 = 0(\Imz\°°) 

(3.7) hD z U{z) = A(z)U(z) + L 2^ L 2(\ lmz\°°) . 

Proof. To see ( |3.5| ) we write 

h^-\\U{t + is)v\\ 2 = 2Re(U(t + is)A(t + is)v, U(t + is)v)) < C\\U(t + is)\\ 2 \\v\\ 2 , s > . 
ds 

Let us now take v with ||u|| = 1 so that, by integration, 

\\U(t + is)v\\ 2 < \\U(t)\\ 2 + % f \\U{t + ia)\\ 2 do- . 

h Jo 

Since this holds for every v with ||u|| = 1 we can replace the left hand side of the inequality by 
\\U(t + is)|| 2 , and the standard Gronwall inequality argument shows that 

\\U{t + is)\\ 2 < Ce Cs/h , 
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which is the desired bound. Putting V(t + is) = d z U(t + is) we have 

hd s V(t + is) = U(t + is)d z A(t + is) + V(t + is)A(t + is) = V(t + is)A(t + is) + 0(s°°) , 

|s| < Mogfr- L , V(t + is)\ s=0 =0, 

where the initial conditi on c ame from the e qua tion o n th e real axis: hD t U(t + is)\ s= o— U(t)A(t). 
As in the argument for (3J5), this implies (|3.6|) and (3/7). □ 

Our definitions of pseudo-differential operators and of (the special class of) fo-Fourier Integral 
Operators were global. It is useful and natural to consider the operators and their properties 
microlocally. We consider classes of tempered operators: 

T : C°°(X) -*-C°°(X), 

and for any semi-norms || • ||i and || • H2 on C°°(X) there exists Mq such that 

|[T U || 1 = O(/i- M0 )|H| 2 . 

For open sets, V C T*X, U C T*X, the operators defined microlocally near V x U are given by 
equivalence classes of tempered operators given by the relation 

T~T' A(T — T')B = 0(h°°) : V'(X) — > C°°(X) , 

for any A,B e such that 

WF(A)cV, WF(B)cU, 

V<sV<£T*X, U<gU$=T*X, U,V open. 

The equivalence class T, ft,-Fourier Integral Operator associated to a local canonical graph C if, 
again for any A and B above 

ATB £ I°(X x X;C'), 

where C needs to be defined only near U x V. 

We say that P — Q microlocally near U X V if APB — AQB = 0^2^ L 2 {h°°), where because 
of the assumed pre-compactness of U and V the L 2 norms can be replaced by any other norms. 
For operator identities this will be the meaning of equality of operators in this paper, with U, V 
specified (or clear from the context). Similarly we say that B = T _1 microlocally near V X V, if 
BT = I microlocally near U xU, and TB = I microlocally near V X U. More generally, we could 
say that P — Q microlocally on W C T*X x T*X (or, say, P is microlocally defined there), if 
for any U,V, U x V C W, P — Q microlocally in U x V. We should stress that "microlocally" 
is always meant in this semi-classical sense in our paper. 

In this terminology we have a characterization of local /i-Fourier Integral operators, which is 
essentially the converse of Egorov's theorem: 

Lemma 3.4. Suppose that U = 0(1) : L 2 (X) -> L 2 (X), and that for every A G *°'°(A") we 
have 

AU = UB, Be$J°(I), a(B) = K*a(A) , 

microlocally near (mo,ran) where k : T*X — > T*X is a symplectomorphism, defined locally near 
too. K ( m o) — m o- Then 

U E Ifl(X x X; C') , microlocally near (too, Too), 

C = {(x,€;y,-rj) ■ (y,rj) = «(»,0}- 
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Proof. From Lemma 3.2 we know that there exists a family of local symplectomorphisms, n tl 
satisfying K t (rrio) — too, and n\ = k, kq = id. Since we are working locally, there exists a function 
a{t), such that n t is generated by its Hamilton vectorfield -ff a (t)- Let us now consider 

hD t U(t) = U(t)A(t), U(1) = U, 0<<<1. 

The same arguments as the one used in the proof of ( |3.3| ) shows that U(0) satisfies 

(3.9) [U(Q), A] = O(h) , for any A 6 tff(X). 

In fact, we take V(t) with V(0) = Id microlocally near (mo, mo), so that 

AU(t)V(t) = Ui^Vit^V^BVit)) = U(t)V(t)A + O(h) , 

where we used Egorov's theorem and the assumption that cr(B) = K*a(A). Putting t = gives 
(3.9). By Beals's Lemma ||, Prop. 8. 3] we conclude that U(0) £ 5'°' (X), and hence U is a 
microlocally defined /i-Fourier Integral Operator associated to k. □ 



If the open sets U or V in (3^) are small enough, so that they can be identified with neighbour- 
hoods of points in T*R n , we can use that identification to state that T is microlocally defined 
near, say, (m, (0,0)), to £ T*X, (0,0) £ T*R™. An example useful here is given in the next 
proposition. 

By Darboux's theorem we know that if p is a function with a non- vanishing differential then 
there exists a local canonical transfomation k such that n*p — £i where £i is part of a coordinate 
system in which the symplectic form is the canonical one d({£,,dx)). The quantum version is 
given in 

Proposition 3.5. Suppose that P £ ^>^ k (X) is a semi-classical real principal type operator: 
p = cr(P) is real, independent of h, and 

p = 0=> dp^O. 

For any m a £ E p ^ f {m £ T* X : p(m) = 0} there exists a local canonical transformation 
k : T* X — > T*R" defined near ((0,0), Too), and an h-Fourier Integral Operator, T, associated to 
its graph, such that 

«*6 = p, 

TP = hD xl T , microlocally near ((0, 0), to ) 
T -1 exists microlocally near (too, (0,0)) . 

For the reader's convenience we outline a self-contained proof of this semi-classical analogue 
of the standard C°° result g, Proposition 26.1.3']. 

Proof. By assumption dp(m a ) ^ 0, and consequently Darboux's theorem gives k with the desired 
properties. Lemma |3.2| then gives us a family of symplectic transfomations K t . If To = U(l), 
where U(l) was defined using the family K t , then (3.3) shows that T P — hD xl = E £ vp -1 * 
microlocally near (0,0). Hence we look for A such that hD Xl + E = AhD Xl A^ 1 , microlocally 
near (0,0). That is the same as solving 

[hD Xl ,A] + EA = , 

Since the principal symbol of P is independent of h, same is true for the principal symbol of E, 
e. Hence we can find a £ S°<°(T*M. n ), independent of h, a(0,0) ^ 0, and such that 

-{£i> a } + ea = 

i 
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near (0, 0). Choosing Ag with the principal symbol a we can now find Aj G ^Z 3 ' (T*M") so that 

[hD Xl , A + A 1 + ■ ■ ■ A N ] + E(Aq + A l + ■ ■ ■ A N ) e ^ N '°(T*R n ) . 

We then put A - A x +A 2 H hAjyH which is elliptic near (0, 0), and finally T = v4 _1 T . □ 

Using the proposition we can transplant objects related to P to the much easier to study objects 
related to hD Xl . In particular, we can microlocally define 

(3.10) ker mo (P) d = T^ker (hD xl )) , ker (hD xl )) = {u G P'(M") : hD xl u = 0} . 

Since ker (hD Xl )) can be identified with V(M n ~ 1 ) we can also identify ker mo (P) with P'(R n_1 ), 
microlocally near (too, (0, 0)): 

(3.11) K : V'iW 1 - 1 ) -^ker mo (P), K = T~ V , tt : x ^ {x 2l • ■ ■ x n ) ■ 

4. Quantum time and quantum monodromy 

Let P(z) G C°°(/ 2 ; \l/°' m (X)), / = (—a, a) C M, be a smooth family of real principal type 
operators, with principal symbols, p(z), independent of h. We will assume that 

S p(z ) = f {m G T*X : p(z)(m) = 0} <s T*X , for z £ I 
P(z) is formally self-adjoint for z € I. 

We assume that mo(z) is a smooth family of periodic points of H p ^, with the minimal periods 
T{z) also smooth in z, and the orbits 7(2): 

exp(T(z)7J p ( z ))(TOo(2:)) = ma(z) , 7 = f {exp(tH p ( z ))(m (z)) : < t < T(z)} . 

When no confusion is likely to arise we may drop the dependence on z in the notation. 
Let flbea neighbourhood of 7(0) in T*X, 

ft ~ 7 (0) xB R2 „-i(0,e), 

and we assume that for z G /, the orbits 7(2:) are also contained in ft. We now introduce a 
covering space of this tubular neighbourhood 

ft ~ K x B R2 n-i(0,e) , tt : ft — >ft, 

with the lift of p(z) denoted by p{z), and we will use the same notation for other objects. 
We start with the following 

Lemma 4.1. The tubular neighbourhood, ft, 0/7(0), can be chosen small enough, so that the 
cover ft contains no closed orbits of Hp( z )y z G [—5,(5] C /, for some small S > 0. 

Proof. Let m 1— > t(m) be a smooth function on ft with the property that £(exp(tiJp( ))) = t, and 
that d=Tr*dt, where dt is a well defined one form in ft. Then Hp^t > on the lift of 7, and 
by shrinking ft if necessary we conclude that Hp^t > on ft. By the periodicity and and a 
compactness argument we conclude that this holds for replaced by z £ [—6,5]. Hence there are 
no closed orbits of Plp( z ) in ft. □ 

We will now replace ft by a finite part: ft ~ [— L, L] x B K 2™-i (0, e), L ^> T. 
A classical time function, q{z) G C°°(ft;R), on ft is defined as a solution of 

(4.1) d z p(z) = -{p(z),q(z)}. 



In view of Lemma 4.1 this equation can be solved (strictly speaking that may involve shrinking ft 
further depending on the initial data, but for simplicity of exposition we will ignore this point), 
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and we can in particular consider solutions satisfying q(rho(z)) = 0. In a neighbourhood of 
too = too(0) £ ft we can define q(z) S C°° such that 

q(z) = ir*q(z) , near too, q(z)(mo(z)) = ir : Q, — > O . 

We clearly have <9 z p = {p, q} near mo- This defines the local classical time near too. We also 
define the first return classical time near too by demanding that 

q(z) = TT*q (z) , near expT(0)Hp^m a , q o (z)(m (z)) = q{expT{z)H p{z) {m a {z)) . 

An iteration procedure similar to the one recalled in the proof of Proposition |3.5| gives the 
quantum analogues microlocally defined near too: 

(4.2) d z P(z) = - l -[P(z),Q{z)], a(Q(z)) = q(z) , 

(4.3) d z P(z) = -~[P(z),Q (z)], a(Q (z)) = q (z) , 

Replacing Q(z) by (Q(z) + Q(z)*) /2, we can assume that Q(z) is formally self-adjoint. We clearly 
have 

Qo{z)-Q{z) : ker mo P(z) — *• ker mo P(z) . 

Then Q(z) is the quantum time near too, and Qo is t ne first return quantum time near Too- See 
the proof of Lemma 7.4 for further discussion of these objects in the classical context. 

For (z, w) near (0, 0), and microlocally near too, we can solve the following system of equations 

{hD z - Q(z) L )U(z,w) d = hD z U(z,w) - Q(z)U(z,w) = 

(4.4) 

{hD w + Q{w) R )U{z, w) A = hD z U{z, w) + U(z, w)Q(w) = 

with the initial condition £7(0,0) = Id, and with U(z,w) bounded on L 2 (microlocally near 
(too, mo)): the solvability of the system follows from the fact that 

[hD z - Q(z) L , hD w + Q(w) R ] = . 

We easily check that (as always, microlocally) 

(4.5) U(z,z)=Id, U(z,w)U{w,v) = U{z,v), 

and that U(z,w) is unitary. In fact, 

hD z {U{z, z)) = Q(z)U{z, z) - U{z, z)Q(z) = -[U(z, z),Q(z)] , U(Q, 0) = ld, 

and U(z, z) = Id is the unique solution. The other property is derived similarly: 

hD w (U(z, w)U{w, v)) = -U(z, w)Q{w)U(w, v) + U(z, w)Q(w)U(w, v)=Q, 

U(z,w)U(w,v)\ w = z = U(z,v) . 

By varying Too along the orbit of -ff p (o), and by extending Q(z) maximally forward (+) and 
backward (— ), we can define semi-global versions of U{z, w): 

U+(z,w) microlocally on a neighbourhood of the diagonal over 

{eTq>tH p (m ) : -e < t < T(0) - 2e} , 
U-{z, w) microlocally on a neighbourhood of the diagonal over 

{exptH p (m ) : -e < -t < T(0) - 2e} . 

The operators have the following intertwining property: 
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Proposition 4.2. Microlocally near the diagonal over 

{exptH p (m ) : -e < ±t < T(0) - 2e} , 
and for z,w close to 0, we have 

P{z)U ± {z,w) = U ± (z,w)P(w) . 
Proof. We define P^(w) — U(w, z)P(z)U(z,w) and differentiate with respect to w: 

hD w P*(w) = Q{w)U(w,z)P(z)U(z,w) - U{w,z)P(z)U(z,w)Q(w) = -[P i {w),Q(w)] , 

PKw)\ w = z = P{z) , 



that is, P*(u>) satisfies (4.2) and consequently P^(w) — P(w). □ 

By replacing the local quantum time, Q(z), by the first return quantum time, Qo(z) (see 
( fOl) , (|J)), we also define U (z,w), 

Ucj(z,z) — Id, U(j(z,w)Ufj{w,v) — U(j{z,v) microlocally near mo. 

This definition will be useful when we study the quantum monodromy operator. To introduce it, 
we first define the forward and backward propagators: 

I ± (z) : ker mo(z) (P(z)) — V'(X) 

(4.6) I±(z) = Id\ kcima(z) (p( z)) , microlocally near m a (z), 

P(z)I±(z) = , microlocally near {exp(tiJ p ( z ))mo(z) : —e<±t<T(z) — 2e}. 

That the operators I±(z) are microlocally well defined follows from Proposition j3.5|, and "mi- 
crolocally" is meant via the identification of ker mo ( z ) (P(z)) with 2?'(R™ _1 ) as in ( |3.11|) . We fix 

mo = Wo(0) as in the definition of f2 above, and define 

W+ — a neighbourhood of mo in T*X, 

^ ^ = a neighbourhood of exp((T(0)/2)iJ p(0 ))TO in T*X, 

ff-C [J expiiJ p(0) (^+), 

|t+T(0)/2|<e 

noting that for z small enough, we can replace mo, T(0), p(0), by mo{z), T{z), p{z) in this 
definition. This shows that maps ker mo(z) (P(z)) onto ker exp((T(z)/2)ffp(2))(mo(z)) (P(z)), 

microlocally near W- x W + . This means that the left microlocal inverse exists and we can give 
the following 

Definition. The (absolute) quantum monodromy operator 

M{z) : ker mo(z) (P(z)) — ker moW (P(z)) , 
is microlocally defined near W-\- by 

(4.8) I + (z)f = I-(z)M(z)f , f e ker mo ( z )(P(z)) , microlocally near W- . 
The quantum monodromy operator, 

M{z) : 2?'(K™ _1 ) — > 2?'(R rl_1 ) , 
is microlocally defined near (0,0) £ T*R n , &y 

(4.9) M(z) = J ftT(z)" 1 X(z)if(z), 



where K{z) is as in ( 3.11 ). 

The basic properties are given in 
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Proposition 4.3. Let U(z,w) and U&(z,w) be given by (4.2), (4.3), and (4.4). Then the fol- 
lowing diagram commutes (microlocally near mo): 

ker mo( ^(P(u;)) ker ma(w )(P(w)) 



(4.10) \U(z,w) \U (z,w) 

ker„ lo(2) (P(z)) ^ kcr mo(z) (P(z)) 
Choosing K{z), so that K(z) — U{z,w)K{w), we also have 

(4.11) hD z M(z) = [K(z)- 1 (g (z) - Q{z))K{z)]M{z) , 



where we recall that by (4.2) and ( |4.3| ), Qo( z ) ~ Q( z ) '■ ^ ieT m (z)(P( z )) ~* ker r?l0 ( z ) (P(z)) , and 
hence K(z)~ 1 is well defined. 

Proof. We need to show that U&(z, w)M.(w) — A4(z)U(z, w), and since is naturally defined 
using the covering space, we will translate this into a statement there. We can microlocally define 
P(w) on f2 and then, 

I+{w) : ker Ao(to) (P(w)) — s- ker(P(w)) , 

and we define 

.M(w) : ker I ~„ o(tu) (P(u;)) — ► ker oxp(T( , 1 , ) ^ (ra))( „- l0 (, i , )) (P(w)) , 

by restricting I+(w) to a neighbourhood of exp(T(u')iJp( 11 ,))(?fio(ty)). Since for n : Q — ► fi, we 
microlocally have 

7r* : ker cxp(T(u; ) ff!;(ra)) ( t?to(u ,))(P(u;)) — ► ker mo(lll) P(w;) , 

7T*X(iy)7T* = M{w) . 



Using the quantized version of q in (4.1), we also define U{z,w), so that U{z,w)P{w) 
P{z)U{z, w). In particular we have 

U(z,w)I+{w) = J+{z)U{z,w). 

Restricting (microlocally) to a neighbourhood of 

(exp{T(z)Hp {z) )(m (z)),m (z)) E h x O , 

and projecting to Q x fi, we obtain 

U (z,w)M(w) = M{z)U{z,w). 



To sec ( 4.11 ) we first note that differentiation of K(z) — U{z,w)K{w) and the definition of 
U(z, w) gives 

hD z K{z) = Q(z)K(z) . 

We then use the commutative diagram to see that 

K(z)M(z) = U (z, w)M(w)U(w, z)K(z) . 

Differentiating this with respect to z and using the previous equation gives 

K{z)hD z M{z) = (Q (z) - Q{z))K{z)M{z) . 

We then recall that by |h|) and (fO|), Q(z) - Qq(z) : ker lflo(z ) (P(z)) -» ker^ io(z) (P(z)) , and 
hence K(z)~ 1 can be applied to both sides. □ 
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We can define the Poincare map for 7 with primitive period T: 

C : T*R" _1 — > T*R n ^ 1 , defined near (0, 0) , C(0, 0) = (0, 0) , 

as follows: for a neighborhood of mo 6 7, t/o, Uq/ exp(tH p ) can be identified with a neighbour- 
hood of (0, 0) £ T*M n_1 (using the local identification of p with £1, as in the proof of Proposition 



3.5), with [to ] corresponding to (0, 0). The Poincare map is then given by 



C : K-\[m]) .— » K-^IexpfTflpJm]) , 
" "' " [m]e C/o/ exp(iP p ) , k : T*M n_1 -> E7 / cxp{tH p ) . 

It will always be undestood that k chosen here is the symplectic transfomation corresponding to 
K = K{z) in j3ll] ) and 

To study quantum properties of the monodromy operator it is convenient to introduce x £ 
C C °°(T*X) satisfying 

f 1 near {exp(tiJ p(0) (to )) : e < t < T(0)/2 - e} 



(4.13) X \0 near {exp(ii7 p(0) (too)) :e<-i<T(0)/2-e} 

n {to : x( TO ) 7^ 1} n {m : x(m) ^ 0} C W+ U W_ , 

where W± are as in ( |4.7| ), and f2 is a small neighbourhood of 7. If p± = 1 microlocally near W±, 
and p± = near W T , we define 

[P,x]w± = P±[P,x] , 

where we use the same notation for y and Op h (x)- We then have the basic property of the 
quantum flux (sec B ) : 



Lemma 4.4. Let K{z) be in (3.11). Then 



U(z) d ^ K(z)*[(i/h)P(z), X ]w + K(z) : P^R"" 1 ) — > P'^ 1 ) 

is microlocally positive near (0,0) G T*IR" _1 and independent of x with the properties ( 4.13| ). 
//we replace K(z) by K(z)U(z)~z then 

(4.14) K(^)*[(j//i)P(z),x]w+^(-2) = Id microlocally near (0,0) £ T*R™ _1 . 

Proof. We note that if P(z)u = near W+, and x is another function satisfying J4.13j ), then 

X(z)*[(*// l )P(z), X - X]« = ^(*)(x - X> - K(z)*(x - x)P{z)u = , 

since P(z)u = 0, and K(z)* P(z) = (P(z)K(z))* = 0. The positivity also comes from expanding 
the commutator and using Proposition [3.5| : 

(K(z)*[(i/h)hD Xl ,x}w+K(z)u,u) = (d Xl xp+K(z)u, K(z)u) > (pu, pu) , 

where K(z) is the composition of K(z) and T of Proposition |3.5| , and p = 1 in a neighbourhood 
of (0,0) € T*]R™ -1 (we again use the same notation for the function and its quantization). □ 



From now on, our choice of K(z) in (3.11) is made so that (4.14) holds. We only need to check 
that we still have 

K{z) = U(z,w)K(w) . 
In fact, we have in general, in the microlocal sense, 

K(z)*[(i/h)P(z), X ]w + K(z) = 

K(w)*U(w, z)[(i/h)P{z), X ]w + U(z, w)K{w) = 

K(w)*[(i/h)P(w),x\w + K(w) , 
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and the last expression is unchanged if we replace x by X (the quantum flux property used 
before). We also used the unitarity of U{z, w). 

With this choice of K(z) we have the following important and well known 



Proposition 4.5. The monodromy operator, M(z), defined by (4.9) with K{z) satisfying (4.14) 
is microlocally unitary: 

M(z)* = M(z)- 1 microlocally near (0, 0) G T*R n_1 , 

and it is an h-Fourier Integral Operator: 



M(z) e 



■C{z)'). 



where C{z) is the Poincare map ( 4.12 ). 



Proof. We need to show that for v e X>'(K™ _1 ) with WF h (v) in a neighbourhood of (0,0), we 
have 



(4.15) 



(M(z)v,M(z)v) = (v,v) + 0(h°°)\\v[ 



I f we put u = K{z)v, use ( 4.14 ), and the defmtion of M(z), (| 
( 4.15 ) becomes: 



then the left hand side of 



{K{z)*[P{z), X ]w + M(z)u,K{z)- x M{z)u) = ([P(z), X ]w + M(z)u,M(z)u),. 



As in the proof of Lemma 4.4, we see that for < t < T(0)/2 + e, the right hand side of the 
previous expression is equal to (modulo 0(h°°)) 

([(i/h)P(z), (exp(tH p ^ ) )*x}exp{-tH p(0) )w + I-(z)M{z)u,I-(z)M(z)u) , 

which corresponds to moving the su pport of x m the direction opposite to the flow of H p (q) , and 
simultaneously moving W+ so that ( 4.13 ) holds. 

Similarly, for -T(0)/2 - e < t < 0, the right hand side of ( 4.15 ) is equal to 

([(i/h)P{z), (exp(-tH p{0) )*x} CX p(tH pm )w + I+(z)u 7 I+{z)} . 
For t - T(0)/2, exp(±tH m (W+) C W-, and hence 

([(i/h)P(z), (exp(~tH p{0) )*x} cxp (tH p(0) )w + I+(z)u,I+(z)u) = 
([(i/h)P(z), (exp{tH m )*xl M -tH p(0) )w + U-(z)M(z)u,I.(z)M(z)u) , t ~ T(0)/2 , 

from the definition of M(z), ( j4.8| ). But this shows ( 4.15| ) proving the first part of the proposition. 

To see the second part we use use Lemma [3.4] , and the obvious conjugation properties of the 
solution in the model case discussed in Proposition 3.5: going around the closed orbit we obtain 
that the underlying symplectomorphism is given by the Poincare map. □ 

So far we have discussed only the case of z G K. We can now consider almost analytic 
extensions of the operators Q(z), Qoi z ), U±(z,w), I±{z), and M(z). For that we consider a 
complex neighbourhod of I C M: 



h 



h,L 



{z : Rezel, | Imz| < i/ilog(l//i)} . 



The families of pseudo-differential operators P(z), Q{z), and Q?j{ z ) have almost analytic exten- 
sions given by Lemma 3.1, and we use the same notation for them. We then use Lemma |3.3| and 
(4.4) to extend U(z, w), Ucj(z,w), and U±(z,w) to (z,w) S Ih.L x Ih,L- We then have 

P(z)U,(z, w) = U m (z, w)P(w) , (z, w) G I h ,L X Ih,L , 



QUANTUM MONODROMY AND SEMI-CLASSICAL TRACE FORMULA 



17 



microlocally (that is, in particular modulo 0(h°°)). Indeed, for x,w G I x I, and \y\ < 
Lh\og(l/h), we have, as in the proof of Lemma pO , 

d v [P{x + iy)U{x + iy,w) - U(x + iy,w)P(w)] = 0(y°°) , 
[P(x + iy)U(x + iy,w) - U(x + iy)P(w)}] y=0 = . 

Hence we can define 

I±(z) = U±(z,w)I±(w) , (z,w) G I h ,L x I, 

so that P(z)I±(z) = 0. 

To define an almost analytic extension of M(z) we first almost analytically extend the pseudo- 
differential operator K(z)~ l (Q(z) — Qq(z))K(z), and then use (4.11) and Lemma 3.2. In par- 
ticular, Proposition 4.5 gives, 

M(z)- 1 = M(z)* , |Imz| < Lhlog(lfh) . 

5. Grushin problem near a closed trajectory 

As in the previous section we assume that P(z) is self-adjoint for z G M, and denote by the 
same symbol the almost analytic continuation of P(z). Although the inverse of P(z) does not 
normally exist near 7 — 7(0) for all z € I wc will describe P(z)^ 1 in terms of the inverse of a 

microlocal effective Haniiltonian E y(z) = I — M(z). We will do it first for z real and then use 

the extensions of operators U±(z,w) described at the end of the last section to transplant the 
results to complex values of z. 

To do that we follow the now standard Grushin reduction || , and consider the system 

(5.1) V(z)= ( ^Jj^ \ ■ V'{X)xV'{W 1 - 1 ) — » V'(X)xV'{W L - 1 ), 

defined microlocally near 7 x (0, 0), and where the operators R± need to be suitably chosen. 
We will successively build the operator V(z) and its inverse. We start by putting 

(5.2) R+(z) = K(zr[(i/h)P(z), X ]w + , 

and u with Pu = near W+, R+{z)u, is its Cauchy data. Hence u = K(z)v provides a local 
solution to the microlocal Cauchy problem: 



(5.3) 



P(z)u = 
[z)u = v 



To obtain a global Cauchy problem we need to introduce R-(z). To do that we define 

K f {z) = I+(z)K(z) , K b (z) = I-(z)K(z) , 



where the operators I±{z) are defined in (4.6). We recall the definition of the monodromy 
operator: 

(5.4) K f (z) = K b (z)M(z) microlocally near x (0,0). 



We can build a solution of (|5.3|) in fi \ W- by putting 

(5.5) E+(z)v = X K f (z)v + (1 - X )K b (z)v , 
so that in particular, E + (z)v = K(z)v, in W+, and consequently 

(5.6) R + {z)E+(z) = Id microlocally near (0, 0) G T* 
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Applying the operator, and using (p.4[) we obtain 



-P(z)E + (z)v = [(i/h)P(z),x]w-K f (z)v-[(i/h)P(z), X ]w-K b (z)v 
= \{i/h)P{z), x ] w _K b {z)(M(z) - I)v . 
Hence we obtain a globally (near 7) solvable Cauchy problem by putting 

(5.7) 
with 

(5.8) R_{z) = \{i/h)P{z), x]w_K b {z) . 



\P{z)u + R-{z)u- = 
R + (z)u = v 



The problem (5.7) is solved by putting 



(5.9) 



u = E + (z)v, U-=E_ + (z)v, E_ + (z) d =I-M(z), 



(5.10) 



£(z) 



where E + (z) was g iven by (5.5]. 

The definitions (5.2) and (5.S) give V(z) in (5.1). If the microlocal inverse, £(z), exists, it is 
necessarily given by 

' E(z) E+(z) 
E_(z) E_ + (z) 

where E+(z) and E |_(z) have already been constructed. 

It remains to find E(z) , E-(z), and to show that the resulting operator £(z) is the right and 
left microlocal inverse of V(z). For the right inverse, this means solving 

^P{z)u + R-{z)u_ = v 
R + (z)u = v + 

We first introduce the forward and backward fundamental solutions of (i/h)P(z): 

Lf(z) microlocally defined near (fl x £ f2) + , 
L b (z) microlocally defined near (CI X e 0)_ , 

where (fi x e Q)± is given by 



(5.11) 



(nx e n) ± = ( |J {(exp(tH, 

\m6f! 



p (o))m, m) n fi x fl : -e < ±t < T(0) - 2e} . 



To do that we use Proposition |3.5| and the corresponding local forward and backward funda- 
mental solutions: 

/xi 
v(t, x')dt, 
-OO 

L° b v(x) = - v(t,x')dt, 

v e £'(R n ). 

We will now try to build an approximate solution of (i/K)P(z)u — v using L m (z). For that let 
us put 

u = L f (z)(l-x)v. 



Let us also define xt, Xf satisfying ( 4.13 ) and in addition, 

Xb = 1 on supp x n W + , x = 1 on supp Xf H W + . 
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We now put 

u = L f (z)(l - X )v , 

where we can think of u as being microlocally defined on the covering space of f2, 51 (see the 
proof of Proposition |4.3|). Hence, P{z)u = to the right of the support of 1 — x ( m the direction 
of the flow), in particular on the support of Xf ■ Hence, to the right of the support of 1 — Xi 

u = K(z)K(z)*[(i/h)P(z), X f]w + u 
= K(z)K(z)*[(i/h)P(z), X f]w + Lf(z)(l - X )v ■ 



If we use the notation from the proof of Proposition 4.3 and put Kf(z) = I + (z)K(z), then in 
the forward direction of propagation past the support of 1 — we have in CI, 

(5.12) u = K f (z)K(zy[(i/h)P(z), x f ]w + Lf{z){\ - x )v . 

Similarly, if u = L b (z)xv then, left to the support of x, we have P{z)ii = 0, and we can extend 
u father left, microlocally in Q: 

(5.13) u = K b (z)K(z)*[(i/h)P(z),Xb\w + L b (z) X v . 

We can think of u and u as multivalued in Q and we will define, near W-, 

Lf fV = second branch of u near W- , 
Lf,bV — second branch of u near W— ■ 

With this notation we put 
(5.14) 



u =E (z)vt { \ L T h[ (\ XV t^ { %-} )V , 
\ L b {z)xv + (1 - x)Lbb{z)X 



outside W- 



v) X v + Lf(z)(l-x)v + X L ff( 1 -X) v inW- 

An application of (i/h)P(z) gives 

(i/h)P(z)u =v- [(i/h)P(z),x}w-L hb (z) X v + [(i/h)P(z), X ] w _L ff (z)(l - X )v , 

and using (fyj) and ( ^13|) (where we now drop the ~ as we are taking the second branch of 
and u, and the definition of M(z), we get 

^P(z)u = v-[{i/h)P{z),x]w.K b (z)K(z)*[(i/h)P(z), X b]w + L b (z)xv 

+ [(i/h)P(z), X }w.K f (z)K(zr[(i/h)P(z),Xf]w + L f (z)(l-x)v 
= v - [(i/h)P(z), X ]w.K b (z) (K(zyi(i/h)P(z), X b}w + L b (z)xv 
- M(z)K{zy[(i/h)P(z),Xf}w + L f (z)(l-x)v) • 

In other terms, 

(5.15) ^P{z)E Q (z)v + R-(z)E ^(z)v =v, 
where we defined E (z) by (|]L4|) and 

(5.16) Eo,-(z) = K(z)*[(i/h)P(z),Xb}L b (z)x- M(z)K(zy[(i/h)P(z) lX f}w + Lf(z)(l - X )v . 
If we now put 

E(z) d = E (z) - E + (z)R + (z)E (z) , E_(z) d = Eq,-(z)v - E_ + (z)R+{z)E (z) , 
then £{z) given by ( 5.10| ) is a right microlocal inverse of V(z). 
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To show that it is also a left inverse, we observe that 

V{z)* = y ~R h ( z y^ R q^* ) : &{X) x V'iR"- 1 ) — > V'{X) x X>'(K™ _1 ) , 

is microlocally dehned in the same region as V{z) and is essentially of the same form but with 
W + replaced by W- and x by 1 — X- 

R + (z)* = [(i/h)P(z), X }w + K(z), 
R-(z)* = K b (z)*[(i/h)P(z), x ]w- ■ 

To see this we first note that 

K h {z)*[{i/h)P{z), X ]w-K b {z) = -Id. 



In fact, as in the proof of Proposition 4.5, (4.14) is invariant under the change of x an d W±, as 
long as (4.13) hold. In particular, for < t < T(0) — e, 

K b (z)*[(i/h)P(z), (exp(tH p{0) )*x} CX p(-tH p(0) )w + K b (z) = Id. 

For t ~ T(0)/2, W + is moved to W-, and (exp(— tH p ^)*x satisfies the properties of 1 — X- 
Hence, using the idependence of Xi 

K b {z)*[{i/h)P{z), (1 - x)]w-K b {z) = Id v>{Mn -i) , microlocally near (0,0) £ T*M" _1 . 

If we now replace K(z) by K b (z), then K(z) plays the role of K b (z), and this proves that R+(z)* 
is the same as — R- (z) with W+ and W- switched and x replaced by 1 — x- 

Hence, a similar argument to the one used for the construction of £(z) shows that Viz)* has 
a right inverse, 

F{z) F + (z) 



Then T(z) is a left inverse of V(z), and the usual argument {J-{z) — J 7 (z)V(z)£(z) = £(z), 
microlocally) shows that it is equal to our right inverse. 

Remark. By constructing part of the left inverse directly we can arrive at a simpler expression 
for E-(z): 

(5.17) = -(M(z)K f (z)* X + K b (z)*(l - x )) , 

and it is useful to have it. To obtain it we will directly solve the problem 



(5.18) 



E-(z)(i/h)P(z) + E-+(z)R+(z) = 
E-{z)R-{z) = /dp'(r-i) 



Motivated by the structure of E + (z) and the fact that R-(z) is close to being an adjoint of R+(z) 
(if it were, then E—(z) would simply be the ajoint of E+(z)), we put 

(5.19) E_(z) = -{M{z)K f {z)* X + K b (z)*(l - x )) ■ 

We now compute 

-E_{z)R^{z) = (M(z)K f (z)*x + K b (z)*(l - X ))[(i/h)P(z),x]w-K b (z) 
= K b (zy[(i/h)P(z),x}w-K b (z). 

To analyze the last expression, we note that K(z), in the definition of K,(z) was chosen, in 
Lemma 4.4, so that K(z)*[(i/h)P(z), x]w + K{z) = Id. As in the proof of Proposition |4.5| , this is 
invariant under the change of x an d W±, as long as (4.13) hold: for < t < T(0) — e, 

K b (z)*[{i/h)P(z), (exp(tH p(0) )*x\cxp(-tH p(0) )w + K b (z) = Id. 
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For t ~ T(0)/2, W+ is moved to W-, and (exp(— tH p r ))*x satisfies the properties of 1 — %. 
Hence, using the independence of y, 

K b {z)*[{i/h)P{z), (1 - X )]w-K b {z) = Id v , (Rn -i) , microlocally near (0,0) e T*K" _1 . 

This shows that _E_(z)i?_(z) = Id and we need to verify the first identity in ( |5.1g| ). For that we 
use K.(z)*P(z) = 0, M(z)K f (z)* = K b (z)*, near (0,0) x VF_ C r*M" _1 x T*Jf, to obtain 

-E_(z)(i/h)P(z) = (M(z)K f (z)* X + K b (z)*(l - X ))(i/h)P(z) 

= M(z)K f (z)*[ X , {i/h)P{z)] w+ - K b (z)*[ X , {i/h)P{z)] w+ 
= K{z)*[{i/h)P{z),x]w + -M{z)K{zT[{i/h)P{z),x]w + 
= (l-M(z))R+(z)=E- + (z)R+(z), 



and that establishes ( 5.1 8j ) , so E_(z) = E-(z) and we have ( |5.17|) . 



So far we considered only the case of z e K, and P{z) — P{z)*. Arguing as at the end of 
Sect.||, we see that all the operators occuring in the construction of V(z) and £(z) have almost 
analytic extensions to | Imz| < Lh\og(l/h) for any L. It follows that the extention of £(z) is a 
microlocal inverse of the extension of V(z) modulo | Imzl 00 , which in this neighbourhood of the 
real axis is 0(h°°), that is, it remains a microlocal inverse. The bounds on the continuation of 
£{z) follow from (|3.5[). This gives 



Proposition 5.1. Let P(z) be an almost analytic extension of the self- adjoint family of operators 
P{z) G C°°(/ 2 ; *°< fc p0) ; such that 

The flow of Hp has a closed orbit 7, 
on which p — er(P(0)) = and dp ^ 0. 



Then, there exist operators R±{z), defined in \ Ivaz\ < Lhlog(l/h), such that 

V(z)= ( ) : V'{X) xV'iW"- 1 ) — ► V'(X) x P'(R" _1 ) , 

defined microlocally near 7 x (0,0), has a microlocal inverse there: 



£{z) 



( E%) E E + }%) ) ^ <^ C|Imz|/ ' 1 ) : L\X) x L»(R-i) -» L\X) x ^(R- 1 ) , 

B z V(z) = 0(\lmz\°°) 



and E y(z) = I — M(z), where M(z) is the quantum monodromy operator defined by (4.9). 



Remark. The constant C in the estimate of the norm of £ (z) could be described more explicitely 
if stronger conditions on P(z) were made. If we assumed (6.1) then C could be related to C p in 



(6.11) 



6. Proof of the trace formula 

We can now prove the main result of the paper. We strengthen our assumptions further here 
by demanding that P(z) is a smooth family of operators, self-adjoint for the real values of the 
parameter, and elliptic off the real axis. 
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Theorem 2. Let P(z) G C°°(/ z ; ^^(X)), I = (—a, a) C R, 6e a family of self-adjoint, principal 
type operators, such that T, z — {to : a(P(z)) — 0} C T*X is compact. We assume that 

a{d z P{z)) < -C < , near E z . 
((K " 1^00)1 > , /or|f|>C. 

PFe aZso assume that for z nearO, the Hamilton vector field, H p ( z ), p{z) — a(P(z)), has a simple 
closed orbit j(z) C So with perid T(z), and that 7(2) has a neighbourhood such that 

(6-2) 

meH and exp t H p ( z ) (to) = to , p{m) = , < |t| < T{z)N + e, z G / , to G 7(2) , 

where T(z) is the period of 7(2:), assumed to depend smoothly on z. Let A G \IV (X) oe a 
microlocal cut-off to a sufficiently small neighbourhood 0/7(0). 

TTien if P(z) is an almost analytic extension of P(z), 
almost analytic extension, f G C°°(R), and supp / C (—N(C P — e) + C, N(C P — e) — C) \ {0}, we 

V / /(z//^ z [x(«) d z P(z) P(z)- 1 ] AC(dz) = 



IT 

(6.3) Af-i 



-L V tr / f(z/h)M(z,h) k ^-M(z,h) X (z)dz + 0(h°°), 

Z7rl _ A r_ 1 aZ 



where M(z,h) is the quantum monodromy operator defined in (4.9) wra'f/i i^(z) satisfying (4.14). 
T/ie constant C p > 0, m £/ie condition on f depends on p{z) only and is given in (6.11). 



We observe that the left hand side of (6.3) is independent of the choice of the almost analytic 
extension of x : if x is another extension then, then x — X — 0{\ Ithz\°°). In view of Lemma 
6.1 below, 

(X(z) ~ xHz)) d z P(z) P(z)- 1 

is smooth in z, and 0{\ Imz|°°) . By Green's formula and holomorphy of /, the corresponding 
integral vanishes. 

As described in Sect.|| Theorem [j] is an immediate consequence of Theorem ^[ 
Before proceeding with a proof we remark that we can assume that 

P(z)e*°'°(X), 

since P(z) can be multiplied by an z-independent elliptic B G flffi k (X), without changing (6.3). 
We start with a lemma which justifies taking the traces in ( |6.3| ): 

Lemma 6.1. Under the assumptions of Theorem^, P{z)~ 1 exists in U\I, where U is a complex 
neighbourhood of J (s /, and 

ll-PW" 1 !! < Cllmzl" 1 , < \lmz\ < 1/C . 

Proof. Let ip — ip w (x,hD; z) be a microlocal cut-off to a a small neighbourhood of E z . Let us 
put v = P(z)u, so that (semi-classical) elliptic regularity gives 

(6.4) ||(1 -^)«|| < C\\v\\ + Oih^M . 

For complex values of z we write 

P(z) = P(Re z) + Im z Q(z) , 
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where P(Rez) is self-adjoint and a(Q(z)) > 1/C > near S z . This shows that 
(6.5) Im{P{z)ipu,il>u) =Imz Re(QO)^«, ipu) >Imz (||V>u|| 2 /C - OQi°°)\\u\\ 2 ) , 

where we used the semi-classical Garding inequality. 
We also write 

Im(P(z)u, u) — lm(P(z)ipu, ipu) = Im z ({Q(z)u, u) — (Q(z)ipu, ipu)) = 
Imz 0(1)||(1 - ^)u||||u|| - Imz 0(1) {\\v\\\\u\\ + 0{h~)\\u\\ 2 ) , 



where we used elliptic regularity (f^J) in the last estimate. Then, applying ( |6.5| ). 

\\u\\\\v\\ > lm(P(z)ipu,ipu)-Imz 0(1) (\\v\\ \\u\\ + e>(/i°°)|M| 2 ) 
> Imz (Uuf/C - 0(1)||«||||«|| - 0(h°°)\\u\\ 2 ) . 

For small Imz the term ||w|| ||u|| on the left hand side can be absorbed in the right hand side, and 
by adding Imz||(l — ^)m|| 2 to both sides we obtain 

Imz||it|| 2 /C < IMHMI + O (h°° ) Im z\\u\\ 2 , 

and that gives 

\H<^-\\v\\, 

Im z 

proving the estimate for P(z) . □ 



Proof of Theorem: Using Proposition 5.1 we can formally write 



P(z)~ L A = E(z)A - E + (z)E- + (z)- L E_(z)A, £L+(z) = I - M{z) , 

microlocally near 17, and for < |Imz < Lh\og{l/ h), with any L. To apply this formal 
expression rigourously, we rewrite the left hand side of (|6.3|) as 



(6.6) 



itr / f{z/h)B z [ X {z) dzP{z) P(z)- 1 } A C(dz) 



-Vtr / f(z/h)d z [ X (z)d z P(z)P(z)- 1 }AC(dz). 
* V Jc + 



Then, motivated by the formal Neumann series expansion of (/ — M(z))^ 1 we define 



AT 

def 



(6.7) T+(z) d ^ E(z)A - E + (z) ]T M(z) k E_(z)A , 

fc=0 

so that 

(6.8) P(z)-M = T+(z) + P(z)- 1 R^(z)M(z) N+1 E^(z)A , 

microlocally near 7 and for < Imz| < Lh\og(l/h), for any L. In fact, from V(z)£(z) = Id, 
and E±(z) = I — M(z), we have 

P{z)E + {z) = -R_{z){I - M(z)) , P{z)E{z) = I - R_{z)E_{z) , 
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and hence 

N 



P{z)T+(z) = P{z){E(z)A - E + {z) Y, M(z) k E-{z)A) 



k=0 

N 



= A-R- (z)E- (z)A + R_ (z) {I - M{z)) M(z) k E_ (z)A 

k=0 

= A-R_ (z)M(z) N+1 E_ (z)A , 



which gives ( p.Sp . 

To use this in ( |6.6| ) we need to have the support of the almost analytic extension of the cut-off 
function \ to be contained in the region where | lmz\ < Lh\og(l/h). To do that we follow the 
method of [[?J Sect. 12] by fixing an almost analytic extension of Xi X^j an d then putting 

- - # , / / \ ,( Im« \ ( 1, |t| < 1/2 

x = x L , h = x*i>L, h , M*) = *[ Lhloe{1/h) ) . = ( 0) |t| >i ■ 



By the remark after the statement of Theorem |2|, ( |6.3D is independent of the choice of x an d 
hence we can use xl,k- Since n ow Q{\ Im z\°°) = 0(h°°), the almost analyticity of P(z) also 
shows that the left hand side of (|6.3|) can be rewritten as 



(6.9) itr / f(z/h)5 z x(z) d z P(z) P(z)- 1 AC(dz) 



and this is what we will use from now on. 
We claim that with the choice of x above 



tr / f(z/h)d z x(z) d z P{z) P(z)- 1 A C{dz) = 

(6.10) \ 

tr / f(z/h)d z x(z)d z P(z)T+(z)C(dz) + 0(h L / c ) 



where C is fixed depending on N and supp /. 
To show this we first need the following 

Lemma 6.2. The almost analytic continuation of the monodromy operator satisfies, for z suffi- 
ciently close to 0, and for any L, 

\\M{z)\\ < e -(C P -e)lmz/h + Q( h ao} ^ < Im Z < Lh\og{l/h) , 

(6n) || JV^T(^)— 1 1| < e < c »- e > In, */ h + 0(/i 00 ) ) -LMog(l//i) <lmz<0, 

,T(0) 

C v = - / a(d z P(z))(exp(tH pi0 ))(m )dt, 
Jo 

where e > can be taken arbitrarily small by shrinking the neighbourhood of"/. The constant C p 
is positive thanks to (|5.l|) . 



Proof. We use the differential equation (4.11) and observe that for z real, and mo{z) £ 7(2), 
a(X(z)- 1 (g (z) - Q(z))K(z))(0,0) = - / a(d z P(z))(exp(tH Az) )(mo(z))dt . 
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Hence, writing z = x + iy, < y < Lh\og(l/h), and B(z) = K{z) 1 {Q(z) — Q(j{z))K{z), we 
have, for v £ D'(M n_1 ), with WF(v) close to (0,0), 

h— (\\M(z)v\\ 2 ) = h4-(M(z)v,M(z)v) 
ay ay 

= ih{d z - d z )(M(z)v, M(z)v) 

= -{B(z)M(z)v,M(z)v) - (M(z)v,B(z)M(z)v) + O(\lmz\ oc )\\v\\ 2 
= -({B(z) + B(z)*)M(z)v, M(z)v) + 0(| ImzDIMI 2 ■ 

The Carding inequality now shows that for x small enough, 



h- (\\M(z)vf) < -(C p - e)\\M(z)v\\ 2 + 0(y°°)\\v\\* • 



Since by Proposition |j| ||M(a;)i>|| 2 = ||u|| 2 (l + 0{h°°)), the lemma follows. 
Proof of ( |S.1C| ) ; By ( |6.6| ) and (^) we need to estimate 

(6.12) tr / f(z/h)d z XL,h^)9 z P(z)P(z)- 1 R_(z)M(z) N+1 E_(z)AC(dz), 



□ 



where by Lemmas 6.1 and p.2| we have 

\\M(z) N+1 \\ < \\M(z)\\ N+1 < e ~( c p~^ N + 1 ) 1J Tr , < Imz < Lhlog(l/h) , 
IIP^HI/Ilmzl. 

All the operators coming from V(z) and S(z) are bounded by exp(C| Imz\/h), and if supp/ C 
[-b + C,b-C], then 

\f(z/h)\ <Ce^ c ^. 
Using the definition of XL,h, the above estimates, and the characteristic function 

PL,h(t) — %/ilog(l//i)/2<i<iftlog(l/ft) j 

we can bound ( |£U2| ) by a constant times 

/ (I<9 Z X # (2)| + ( J L/ilog(l/^))- 1 i OL,/ i (Imz)) |Imz|- 1 e^( b -( c "- e )( JV + 1 »/:(dz) < 

C/T n / |Imz| M e^( b -( c ^ e )^ +1 »/:(dz) 

JO<Iniz<Lh\og(l/h) 

+ Ch~ n {Lh\og{l/h)y 2 h« c ^ N+1 ^ L < C' N h L ' c - n , 

where C > is fixed. □ 

With ( 6.1 0| ) established, we have to study the leading term on its right hand side which we 
rewrite using the definition (6.7) and the cyclicity of the trace: 



-tr / 5 z (xL,h)d z P(z) E(z) A f(z/h)C(dz) 



(6.13) 



-tr / d z (xL,h) \ y j M{z) k E^(z)Ad z P{z)E + {z)\f{z/h)C{dz 



2(1 



J. SJOSTRAND AND M. ZWORSKI 



Since all the operators are almost analytic (in particular d z E,(z) = 0(/i°°) on the support of 
Xlji) and f(z/h) is holomorphic, we can apply Green's formula and reduce the first integral to 
an integral over the real axis: 

-^tr [ X (z)f(z/h)d z P(z)E(z)Adz. 
2m J R 

To analyze the sec ond term ( with integration still over C+) we see that the explicit expressions 
E-(z) and E + (z), and (|5.17D , show that 

E_(z)Ad z P(z)E+(z) = Ai(z)M(z) + A 2 (z) , A 3 {z) e C°°(/ z ; (R"" 1 )) . 

To analyze the contributions to the trace we need the following simple 

Lemma 6.3. Suppose that A E V%-°°(R n - 1 ) and U G /"(M™ _1 x R n_1 ; C), satisfy 

m e WF h (A) => {m,m) i C . 

Then, for X e C£°(M. n ~ l ), 

ti X AU = 0(h°°). 

Proof. It is clear that if xi,X2 £ C C °°(T*R"- 1 ) have disjoint supports then tr X iAUx2 = 0(h°°), 
where we denoted the quantizations by the same symbols. Using the hypothesis, we can write 
X-AU as a sum of negligible terms (0(h°°)), and of terms of that form. □ 



The assumption ( |6 . 2| ) implies that the main contribution (modulo 0{h°°) as usual) to the 
trace of 

M{z) k E^{z)Ad z P{z)E+{z), k>0, 

comes from an arbitrarily small neighbourhood of the fixed point, (0, 0) G T*R n_1 , of C(z), the 
canonical relation of M(z). We can therefore replace A by 1 and introduce a microfocal cut-off, 
p w , to a neighbourhood of (0, 0): 

tr M(z) k E_(z)AdP(z)E+(z) = M(z) k E^(z)dP{z)E+(z)p w + 0{h°°) , k > 0. 

For k — the same discussion is valid for the contribution of M(z)A\(z) in E-(z)dP(z)E + (z), 
but for the pseudo-differential contribution, A 2 (z), we need to use the support assumption on /: 
^ supp /. We write 



tr / f(z/h)d z x(z)A 2 (z)p w £(dz) = tr / f(z/h) X (z)A 2 (z)p w dz = 0{h°° 

J C+ JR 

by the standard argument: put g(z) = tr A 2 (z)p w , so that by PlancherePs theorem 

(6 ' 14) = /* M+1 / f(hC)/(hC) M C M 9(-0d( 

Jn 

= 0(h M+1 ). 
Hence the second term in ( 6.1 3| ) becomes 

(6.15) --tr / B zX [ y j M{z) k E^z)d z P(z)E + {z) j p w f(z/h)C(dz) , 

where p w is a microfocal cut-off to a neighbourhood of (0,0). 
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We recall that when R± are independent of z, the following standard formula holds: 
d z M(z) = -d z E-+(z) = E-(z) d z P{z) E+(z) , 
as is easily seen from d z £ = —£d z V£. In the general case, the same argument gives 
d z M(z)=-d z E_ + (z) 

((K -'" =E-(z) d z P(z) E + (z) + E-+(z) d z R+{z) E+(z) + E.(z) O z R-(z) E-+(z) . 



Inserting this we obtain the following expression for (6.15): 

1 r N 

- - tr / f(z/h)dx(z) V M(z) k d z M(z)p w C{dz) 

7T /it Z — ' 

• yL + fe=0 

AT 

(6.i7) tt y c+ ^ o 

JV 

-tr / /(V/i)^)X) M (*) fc£ -(*) fl * B -(*)( 1 - M (*)) _ V u 'Af(*) fc A d2! ) 



77 "W- fc=0 
= Ji + J 2 + .h 
By Green's formula 



Ji=~f]tt / f(z/h) X (z)M{z) h d z M{z)p v 'da, 
2717 fcO Jr 



which is a term appearing in ( |6.3| ). We want to show that the remaining two terms, J2, J3, are 
negligible. 

To see this we need 

Lemma 6.4. We have 

d z R + (z)E + (z) , E_(z)0,J2_(z) e C°°(7 Z ; . 



Proof. From the definitions (5.2), (5.5), and from (4.14) we see that 

d z R+(z)E + (z) = d z (K*(z)[(i/h)P(z), X ]w + ) K{z) = -K*(z)[(i/h)P(z), X ]w + d z K(z) . 
From the proof of Proposition [Ij^ we recall that 

hD z K{z) = ~Q{z)K{z) , 

and hence 

d z R+(z)E + (z) = (i/h)K*(z)[(i/h)P(z),x] W+ Q(z)K(z). 

This expression is microfocal near (0,0) and as far as K (z) is concerned microfocal near (mo, (0,0)). 
Hence we can use a model given in Proposition 3.5: P(z) = hD xl (the microfocal ^-dependent 
conjugation will not affect the uniform pseudo-differential behaviour), and 

K(z)u(x 1 ,x') = -J— 1 J e^^-*'^^ h a z ^,tf)u{^H, 

where we used local representation of the h- Fourier Integral Operators (see the proof of Propo- 
sition 7J3 below for the derivation of a local representation). After composing the operators, 



and applying the stationary phase method we arrive at the following expression for the kernel of 

d z R + (z)E + (z): 

1 



{2nh) n - 1 
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which by a standard "Kuranishi trick" argument (see the appendix) shows that we get a smooth 
^-dependent family of pseudo-differential operators. □ 



In J 2 we can replace Y^k=o M(z) k (l - M{z)) by 1 - M(z) N+1 . As in the proof of ( |6.10| ), we 
show that the term corresponding to M(z) N+1 is negligible. The remaining term is transformed 
to an integral over R: 

f f(z/h) X (z)tr(d z R + (z)E + (z)p w )dx, O^supp/ 



which is negligible by Lemma 6.4 and ( |6.14 ). Similar arguments then apply to J; 



To summarize, we have shown that 

tr / f(z/h)8,x(z) d z P(z) P(z)- 1 AC(dz) = 

JC+ 

ic 1 N r 

— tr / X (z)f(*/h)d,P(z)E(z)Adz-— Y> / f(z/h) X (z)M(z) k d z M(z)p w dx + 0(h°°) , 

zl7rl /l7rl fc=Q JR 



We proceed in a similar way for the integral over C- in (6.6). We write I — M{z) = —M (z)(I 



M(z) ), and motivated by the resulting formal Neumann series put 

N 

T N {z) ^ E(z)A + E + (z)J2(M(zy 1 ) k E^z)A . 

fe=l 

The same arguments apply and Green's formula gives 

i r i r N+1 

-tr x{z)f{z/h)E{z)Adz-—M X(z)f(z/h)Y] M(z)- k d z M(z)dz + 0(h°°) 

Jib ^Ttl ./it k—\ 



When we now add the contributions from the integrations over C± we see that the integrals 
involving E{z) cancel and the remaning terms give d6.3| ) 

□ 

7. Trace formula for non-degenerate closed trajectories 
We say that a closed trajectory y(z) of P(z) N-fold non-degenerate if 

(7.1) det(I-(dC(z) ma(z) ) k )^Q, 0^\k\<N, 

where C(z) is the Poincare map of j(z), ( j4.12| ). When this holds our theorem translates 
into the standard semi-classical trace formula, generalizing the Gutzwiller, Balian-Bloch, and 
Duistermaat-Guillemin trace formulae. 

We start by a general discussion of traces of Fourier Integral Operators. 

Lemma 7.1. Suppose that B, microlocally defined near (0,0) £ T*R™ is given by 

(7.2) Bu[x) = -1— J f e m*«)-mVh b ( x , m h)u(y)dydn 

where (f>(x,rf) is defined near (0,0), a is a classical symbol of order 0, supported near (0,0), 
</>'(0,0) = 0, and 0^(0,0) ^ 0. The corresponding canonical transformation is given by 

(7.3) K:(^ n (x ) Tj),ri)^(x,(f/ x (x,Tf)). 
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It is defined between two neighbourhoods of (0, 0) and we assume that (0, 0) is its only fixed point 
there, and 

(7.4) det(d/s(0,0)-l)^0. 
Under these assumptions 

(7.5) tr£ = 22 - — LJ — £ L^i s = sgn ™ , T JP 

|det^det(dK(0,0)-l)|i l<x-l 

where the signature of a symmetric matrix A, sgn A, is the difference between the number of 
positive and negative eigenvalues. 



Proof. The fact that (7.3) defines a smooth map is equivalent to the assumption that 

(7.6) detC ; ^0. 

Here and in the following, second derivatives of (f> are computed at (0,0) if nothing else is specified. 
The differential, dn(0, 0), is the map (Sy,S„) <— * (6 X ,S^), where 

S £ = <t>l v 5 v + CA ■ 

Here we can express S x and 6^ in terms of S y , 5 V and it follows that dn(0, 0) is given by the matrix: 

(7.7) d/t(0,0) = 



w 1 -w-v; 



" - 6" (6" 



We find the following factorization: 
d/c(0,0) - 1 = 

( 7 - 8 ) ( -CO" 1 oWo 1W1 C-i 
o i M i o \ o i M - 1 



TjX >VV ' 



In particular, 

(7.9) det(d«(0, 0) - 1) = — ^ det 



b" 6" - 1 

XX YXT\ 



, y> vx 1 § m 



Here 



xx rxr] 



is the Hessian of <j)(x,ij) — xr\. The stationary phase method applied to the trace of (7.2) gives 

tiB = (H( &-i c" 1 ))"" , « M)+ow ) eW *' 

Here we choose the branch of the square root of the determinant on the set of non-degenerate 
symmetric matrices with non- negative real part which is equal to 1 for the identity. Using (7.9), 



we get (7.6). □ 



To give a geometric meaning to the signature s appearing in (7.5) in terms of a Maslov index we 
first recall the definition of the Hormander-Kashiwara index of a Lagrangian triple: let Ai, A2, A3 
be Lagrangian planes in a symplectic vector space (V,lo), and put 

(7.10) s(Ai, A 2 , A 3 ) = sgn Q(X U A 2 , A3) , 

where Q(Ai, A 2 , A3) is a quadratic form on Ai © A 2 © A3 given by 

Q(Ai, A 2 , A 3 )(ui © v 2 © v 3 ) = u>{v 1 ,v 2 ) + uj(v 2 ,v 3 ) + w(t> 3 ,ui) , 
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see Q for a comprehensive introduction. Here we only mention that if Aj's are mutually transver- 
sal, then s(Ai,A2,A3) is the only symplectic invariant of such three Lagrangian planes. It is 
antisymmetric and satisfies the cocycle condition. 
We then have 

Lemma 7.2. Let V = T*W l x T*R™ with the symplectic form lj = uj\ — LU2, where lj\ and 0J2 



are the canonical forms on the factors. In the notation of Lemma 7.1, let Tdn be the graph of 
d/c(0,0), A C T*R" x T*R" be the diagonal, and M = {0} © R n © R™ © {0} C T*K™ x T*R™. 
Then 



^xx *Pxr} ^ 

b'L-l <b'L 



s{T dK ,A,M) = -sgn 

Proof. Let us write 

Ao,o)=(^ * 

Then Tdn — {(x, ax + (3rj; (3 l x + 777, 77) : (x, rf) 6 R" x R" }. Since A and M are transversal, || 
Lemma 1.5.4] says that 

S (r dK ,A,Af) = -sgnc;(7r.,.)|r £iK , 

where n : T*R™ x T*R" — > M is the projection along A: 7r(x, £; y, rf) = (0, £ — 77; y — a;, 0). In the 
(x, 77) coordinates on T^, w(7r», •)|r dK is then given by (f2», •), where 

a 13-1 
P-I 7 

which proves the lemma. □ 
As is well known, and as will be seen in the proof of the next proposition, any locally defined 



n = 



Fourier Integral Operator can be represented by (7^). To compute its trace in terms of invariantly 



defined objects we also have to recall the definition of the Maslov index of a curve of linear 
symplectomorphisms - see for more details and references. 

Thus let T(t) C T*R™ x T*R", a < t < b, be a curve of graphs of symplectomorphisms. Choose 
a subdivision a = t Q < t\ < ■ ■ ■ < t^ = b, such that, for all j = 1, • • • , k, there is a Lagrangian 
subspace Mi transversal to T(t) and the diagonal, A, for t € [tj-i,tj]. We now follow Q and 
define the Maslov index of a curve of linear symplectomorphisms as 

(7.11) fx ^ ±£ (s(rfe-i), A, Mi) - s(T( tj , A, Mi))) . 

i=i 

It is independent of the choice of the transversal Lagrangians, Mi, and of the subdivision. 
We can now prove 

Proposition 7.3. Suppose that U(t) is a family of Fourier Int egra l Operators defined using a 
family of pseudodifferential operators, A(t) £ ty°'°(T*X), as in (3.2): 



hD t U(t) + A(t)U(t) = , U(0) = U E *%°(T*X) . 

Let us also assume that at = a{A{t)), the Weyl symbol (with a possible dependence on h in the 
subprincipal symbol part) of A(t), is real and generates a family of canonical transformations: 

—Kt(x,(;) = (Kt)*(H at (x,t)), K {x,£) = {x,£), (x,£)ET*X 
« t (Q,0) = (0,0). 
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If 

(7.12) det(l -<2k t (0,0)) ^ 0, 
then 

(7.13) trU(T) = i^ 1 + °^'-' I ^ M : ,/ \ mm , 

K J |det(d« r (0,0)-l)|i 

where fx(T) is the Maslov index of the curve of linear symplectic transformations dn t (0,0), < 
t < T. 

Proof. Let us first assume that for < t < T 

(7.14) (K t (y,n); (y,rj)) h-> (x(n t (y,r))),x) is surjective near (0,0) . 

We follow the presentation from || Appendix a] . Let at be the Weyl-symbol of A t defined mod- 
ulo 0(h 2 ) (if there is a subprincipal symbol we include it in the principal one and obtain an h 
dependent symbol) . Consequently the influence of the subprincipal symbol will be accounted for 
as an ©(^-dependence in the canonical transformation K(.) Let Kt be the canonical transforma- 
tion generated by H at as described in the statement of the proposition. We can then view n t as 
the canonical transformation associated to U(t) (defined modulo 0(h 2 )) and we claim that 

(7.15) U(t)u(x) = 1 / / e&W'^-v^bit, x, r); h)u(y)dydrj , 



(2irh) n 
where 

(7.16) d t cj>{t, x,rj) + a t {x,d x 4>) = 0, 0(0, x, i]) = x ■ rj , 

so that Kt : ((?,,</>, 77) 1— > (x,d x (f>). The amplitude b has to satisfy 

(hD t + af (x, hD))(eW t ' x >'iV h b(t, x, 77; h)) = 0, 
' " ' ' ' (8t<f> + hD t + e-^ /h af{ Xl hD)e^ /h )(b) = . 

Here the Weyl symbol of e -i<t>/h a w e i<t,/h ig g^^) = a t (x,<f>' x + + 0(h 2 ), and using that 
dt4> = —at(x, d x <t>), we get 

(hD t + Op (ft(x,Z)))b = 0(h 2 ), 
where ft(x, £) — a t (x, 4>' x + £) — a t (x, (f>' x ) (and 77 is just a parameter). This gives 

1 " 

(7.18) (hD t + x Y^^'K <<i>hl), + hD X] o (d^a t )))b = 0, 

1 

for the leading part of b. With v t = ^Zip^,at)d Xi , this can also be written 

(d t + v t + ^divi/ t )6 = 0, 



or 



(7.19) (d t + C Ut ){b (t,x,r ] )(dx 1 A ... A cfe„) 1/2 ) = , 

where C Vt denotes the Lie derivative. If we consider bo(t, x, rj)(dx\ A ... Adxn) 1 / 2 as a half-density 
on n = {(x, dx4>t(x,n))}, then ( 7.19| ) means that 

(7.20) K* t (b Q (t,x 1 n)(dx 1 A...Adx n ) 1 / 2 \ A , ) = (dx 1 A...Adx n ) 1 / 2 \ A . . 



From (7.7) it follows that the restriction of the differential of n t to TA^^ followed by the 
x-space projection is given by 



:S2 
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so ( 7.20| ) says that 
(7.21) 



(det^)V2 



1 . 



We note that det <j% x > for < t < T. 

From QT.16j ) and d(/> t (0,0) = (since «t(0,0) = (0,0)) we see that 

0r(O,O) = -/ ot(0,0)dt. 
Jo 

Applying Lemma 7T and Lemma [7^2] w e obtain ( [7.13 ): under the assumption ( 7.14 ) we only 
need one transversal Lagrangian in (7.11), and we can take M from Lemma 7.2. Then 

/i(T) = ~ (s(T dK0 , A, M) - s (r dKr , A, M)) 

= i( s (A,A,A/)- s (r dKT ,A,M)) 

= 4 S (r dKT ,A,M ) = i S gn(^_ 1 



6" 



I n the case ( 7.14 ) does not hold for < t < T, we have to choose different coordinates in which 
( 7.14 ) holds for tj-i < t < tj. That gives corres ponding Lagrangians M J (defined as M was) and 
the phase shifts add up precisely to give ( |7. 1 1| ) . In fact, we can conjugate U(t) by an h- Fourier 
Integral Operator (so without affecting the trace), so that for t\ — 8 < t < ti + 5 the resulting 
operator is given by 



U(t)u(x) 



(2-Khy 



e ^t,x,r,)- v . v )i^ a . ! m h)u(y)dydr) . 



where we can arrange that b > 0, and that ( [7.12 ) holds with T = t\, t%. We then use Lemma 7.1 
and the geometric discussion above to compute the trace 

r*i 



i - f 1 + Of h))e~ l to 1 a *i (o,o)dt/h 
|det(rf/c tl (0, 0) - 1)| 2 



where 



-sgn 



XX TXT] 

1 



As in Lemma 7.2 



7}X 

we interpret s as 

§ = s(T dKti ,A,M 2 ), 

where Mi was chosen in new coordinates. Comparing it with the previous expression for the 
trace (where we put T = t\ and M = M 1 ), we see that 

v = s(T dKti , A, A/ 1 ) + s(T dKti , A, M 2 ) . 

We now use U (t) to compute the trace at t — ti which in view of the expression for v, and the 
fact that s(r dKo , A, M 1 ) = is 

tvU(t 2 ) = i5(*(ri« ,A 1 M 1 )-s(rj« ll ,A,M 1 )+s(r d « ll ,A,M a )-«(rj« f2 A,m 2 )) x 
(1 + Oih^e'^o 1 a t!(o,o)dt/h 



|det(d« tl (0,0) - 1)| 2 



-a(f/ o )(0,0), 



and comparing with (7.11) see that the power of i is given by the Maslov index for the curve 
dn t (0, 0), < t < t 2 . We can continue in the same way which gives us the final index /i(T). 
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□ 

We now want to evaluate the trace of M(z, h) k M'(z, h), and for this we need to identify the 
Maslov factor and the phase. For this we recall the definition of the classical action: 

(7.22) I(z) d = f £dx . 

The well known relation with the periods is given in 



Lemma 7.4. Let q(z) and qo( z ) be the local time and the first return local time defined in (4.2) 
and (fU). Then 

l-T(z) 

(<lo( z ) - l( z ))U(z) = o-(d z P(z)){exp(tH p{z) )(mo{z))dt, 
Jo 

(goO) -q(z))U( Z ) = j^( z )- 

Proof. The first identity follows directly from the definition and was already used in the proof of 
Lemma |6.2| . 

Since d z p{z) ^ 0, we can write p(z) = c{z){z —p). Hence on p{z) — 0, the equations for q and 
are 

Hpq = -1 , HfAo = _ 1 j 

and consequently {qo( z ) — q( z ))\-y(z) is the period of 7(2), thought of as an orbit of Hp on p = z. 
We now introduce an isotropic submanifold, T, of T*(X x K), where the new variable (on M.) is 
denoted by £ with z its dual variable: 

T = {(m; (C, z)) e T*(M xR) : me 7(2) , C = q( z )(m) , z a < z < z x } . 

The symplectic form dt^Adx + dzAd^ vanishes on T, and hence we obtain from Stokes's theorem: 

I(z\) — I(zo) = / £dx — / £dx 

Jl(zi) J~l(zi) 
pzi rzi 

= I (q(z)(m (z)) - q Q (z)(m (z))) dz = T(z)dz , 

J za J z 

which proves the lemma. □ 

Using this lemma we will be able to identify the phase in the trace of the monodromy operator. 
For that let T(z) be the quantum time appearing in (4.11): 

T(z)^K(z)- 1 (Q(z)-Q (z))K(z), 

so that that formula becomes 

(7.23) hD z M(z) = T(z)M(z) . 



This and Proposition 7.3 show that the phase factor in trM(z, h) k T(z) satisfies 

J'ki z ) = Hlo(z) -?(z))ItM ■ 
In fact, for any family of P(z, h) satisfying the assumptions of Proposition 5.1, we can associate 
to M(z, h) k (not necessarily satisfying the non-degeneracy condition) a phase factor, Jk(z) which 
has to satisfy 

J k (z) = kI{z) + C k , 
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We want to show that Ck = 0. For that we note that if we put P e (z, h) — P(z, h/e) then the 
corresponding Jk(z) is given by Jk(z)e. On the other hand, the action corresponding to P e is 
kl(z)e. Since we can consider P e as another deformation of our operator we must have 

Ve>0, J k (z)e = kl(z)e + C k , 

and for that we need that Ck = 0. 

To obtain the Maslov factor we need to find a family of symplectic transformations interpo- 
lating between the identity and the Poincare map. For that let us fix z and supress dependence 
on z in the subsequent formulae. We want to define a family M{t) : ^'(R™ -1 ) — > 2?'(R™ _1 ), of 
/i-Fourier Integral Operators such that M(0) = Id and M(T) = M. To do this we modify the 
definition of I + in ([0]) to 

I + (t) : ker mo P — > ker exp w p m (P) 
We also generalize the definition of K to 

K{t) : V'(R n ^) — > ker^ ^(P), 



defined using Proposition 3J5 as in ( 3.11 ). We can now define 



(7.24) M(t) = K(t)- 1 1+(t)K(0) : 2?'(M™ _1 ) — > £>'(R™ _1 ) , 

microlocally near (0,0). This family has desired properties and quantizes a curve of local sym- 
plectomorphism re t : 

k 4 = ($ t )-%$ ■ T*R n - x — ► T*R n -\ K t (0,0) = (0,0), 

where &t symplectically identifies a neighbourhood of (0,0) in T*M™ _1 with a submanifold of 
p = 0, St, transversal to 7 at exptH p (mo), and ^ t '■ Sq — > St is the restriction of the flow 
exp(s_ffp) to St- The construction above allows an arbitrary choice of St and $ f . 
We can summarize this discussion in 

Proposition 7.5. Suppose that the orbit 7(2;) is primitive and N-fold non- degenerate in the 
sense that (7.1) holds. Let I(z) be the classical actions defined by ( 7.22| ), and T(z) the period. 



s 



ofl(z). 

Ift, < t < T{z) parametrizes 7(2), let St be a family of submanifolds ofp(z) = 0, transversal 
to 7(2:) at t, $t a symplectic identification of St with a neighbourhood of (0,0) in T*M™ _1 , and 
^t '■ So —> St the restriction of the flow to St- Then for < |fc| < N , 

,m< h \k-Hn Mi h\ e^)e^Wf( g0 (z)~ g (z))| 7(g) /i 

ttM(z,h) hD z M(z,h) = — — — -j (l + 0(h)), 

\(dC(z) mo(z) ) k - 1 1 2 

where v k (z) is the Maslov index of the curve of linear symplectic transformations: 

d^tftSotao) » < t < kT{z) . 

Remark^ The Maslov index v k [z) is a locally constant function of z: it does not change as 
long as (7.1) holds. Its value may depend on the non- unique choices of the identifications <I>t, 



and the transversals 5*4. Since exp(it'/ c 7r/2) is determined uniquely (as it appears in the trace!), 
v is termined only modulo 4. In the case when 7(2;) — > 71(7(2)) is a diffeomorphism, with 
7T : T*X — > X, the natural projection, a choice of transversals submanifolds in the base gives 
natural St's in {p — 0} C T*X. Thus in the case of the geodesic flow on a Riemannian manifold 
v is the index of a closed geodesic. 

The usual semi-classical trace formula for non-degenerate orbits follows from Theorem [l] and 
the following 
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Proposition 7.6. Suppose that the assumptions of Theorem^ are satisfied, M(z) is the quantum 
monodromy defined in (4.9), and in addition the closed orbit 7 = 7(0) is N-fold non- degenerate 
(frTll). Then, for k^O, \k\ < N, and g £ C~(R), we have 

1 + . f^„ m »^ Mfc-i d *,r„ u\..r„\j~- e lkS -< /h+lu -'- k iT 1 g{kT 1 ) 



(7.25) — tr / g(z/h)M(z, hf~ 1 —-M{z, h) x {z)dz = -^^^—B^* + Q{h) , 

J R az I det((aC 7 ) /c - i)| 2 

where T 7 is the primitive period of"/, dC 7 is the linear Poincare map, S-y, the classical action of 
7, and v lt k the Maslov index of k^y. 

Proof. Since P is assumed to be self-adjoint, the subprincipal symbold of P is zero. Let n z be 
the Poincare map and assume that (dK O (0,0)) fe — 1 is non-degenerate. Let T(z) be the quantum 
time appearing in (7.23) above. Using the cyclicity of the trace, we can write the left hand side 
of (|7\25|) as 

(7.26) i-tr / x(z)9(l)M(z) k T(z)f. 

J Ft 

The leading symbol of T(z) at the fixed point is the period T(z) = dl(z)/dz, where I(z) is the 
action along the closed trajectory. Then to leading order, (7.26) becomes 

jH f _ ikl(z)/h i 

(7-27) — / x (z)g(r) r T(«)?, 

V ' 2ttJ M m /i ; |det((d/c z (0,0) fc -l)|3 h' 

where \x is the Maslov index. Write E — z/h, so that 

I M = m +I > m+ o(h). 



Then (7.27) becomes, again to leading order, 

e ifei-(o)/fe T (o) x (o) _ e lkI ^/ h T{0) X (0)g(kT(0)) 



(7.28) ^- J g(E)e ikI 'W E dE- 



det((dK O (0,0)) fc - 1)| 2 |det((ck o (0,0)) fe - l)|i 



□ 



The usual Gutzwiller trace formula for a more general class of operators is given in 



Theorem 3. Suppose that the assumptions of Theorem hold and that in addition 7 is an N- 
fold non-degenerate orbit in the sense that ( f7.l[ ) holds. Then in the notation of Proposition 7.6 
we have 

trf(P/h) X (P)A=±- £ .^ Jy +OW- 
2n k~^ N |det((dC 7 ) fc -7)| 2 



Appendix 

In the classical treatment of pseudo-differential operators, the subprincipal symbol is invariant 
under coordinate changes when the pseudo-differential operators are considered as acting on half- 
densities. This invariance is particularly nice in the Weyl calculus, where the subprincipal symbol 
is contained in the leading symbol - see ^ Sect. 18. 5]. 

For the reader's convenience we present here a self-contained discussion of the analogous result 
in the semiclassical setting. 
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We use the informal notation for sections of the half-density bundles: 

uGC°°(X,r2|) u = u(x)\dx\?, 

a e S°'°(T*X,Ql, x ) a = a(x,£)\dx\i\d£\v , 

which captures the transformation laws under changes of coordinates: 

u{x)\dx\ 5 = 5(x)|oJx| 5 , x = k(x) u(k(x))\k' (x)\ ? = u(x) , 

where for a linear tranformation A we denote its determinant by |^4|. 

We observe that the half-density sections over T*X are identified with functions if we consider 
symplectic changes of variables, and in particular 

(A.l) (x,0^(x,i) = (K(x),( K '(x)n). 



As stated after fl3.1| ) in this paper we considered pseudodifferential operators acting on half- 
densities: 

with distributional kernels given by 

(A.2) K a (x,y)\dxfi\dy\i = J a e^^d^dx^dy^ . 

We will show that 

(A.3) K a (x,y)\dx\$\dy\* = K- a (x,y)\dx\^\dy\i with 5(5, £) = a(x, £) + 0(/i 2 (£>~ 2 ), 



where (x,£) is given by (A.l). 

To establish (A.3) we start with its right hand side using the coordinates x = k(x) and 

1 

(2ith) n 

Making a substition we obtain 



^,tje^-y^/ h di\dx\i\dy\^ 



1 



(2tt/i)™ 7 V 2 
We now apply the "Kuranishi trick" and for that write 

k(x) - n(y) = F(x, y)(x - y) , F(x, y) = k 1 (^t 1 ) + ((x - y) 2 ) , 
(A.4) V J 

K (x) + K(y) = Ki^yj+0((x-y) 2 ). 

We put £ = F(x. y) £ and rewrite the expression above as 
1 f I ' _ ( ( x + y x 



, (K^yY)' 1 ^ + 0(x - yf^j e^-y-V^K^yyr 1 



(2nh) n J V V V 2 

<2£|k'(x)|' |K'(y)| 5 |(ix|2 \dy\* = 

| ^5 (k ((x + y)/2) , (V ((x + y)/2) ( ) ' 1) + 0(x - y) 2 ) ^"^V^ + ?i)/ 2 )P 1 

d^|«'(a:)|^|«'(s/)|a |<ia;|*|dy|= , 

and the terms 0((x — y) 2 ) will influence the symbol only modulo 0((^)~ 2 h 2 ) (by integration by 
parts based on (x — y) exp((a; — y, £) //i) = /i_E>£ exp((x — y, £) //i) ), and hence can be neglected. 
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We now observe that 

W((x + y)/2)\ 2 = \K'(y)\\K'(x)\ + 0((x - yf) , 

and consequently 

Ka{x,y) = 



1 



5 L {(x + y)/2) , (V ((* + y)/2)*) 1 e) + O(h 2 (0- 2 )^ e^-^ldeida^ldi/l* , 



which is the same as K a (x, y)\dx\ s \dy\ % , if 

a(x, o = a (i) , (V (a) 4 ) _1 + o(^ 2 (0" 2 ) • 

This proves (A. 3) completing the appendix. 
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